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Abstract 

Plug conveying along a vertical tube has been investigated through simulation, using 
a discrete element simulation approach for the granulate particles and a pressure 
field approach for the gas. The result is compared with an experiment. The dynamics 
of a plug are described by porosity, velocity and force profiles. Their dependence on 
simulation parameters provides an overall picture of plug conveying. 
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1 Introduction 



A quite common method for the transportation of granular media is pneumatic 
conveying, where grains are driven through pipes by air flow. Practical appli- 
cations for pneumatic conveying can be found in food industry and in civil 
and chemical engineering. One distinguishes two modes of pneumatic con- 
veying: dilute and dense phase conveying. Dilute phase conveying has been 
studied in much detail [1,2,3,4,5,6,7,8] and is well understood. The grains are 
dispersed and dragged individually by the gas flow and the interaction be- 
tween grains is small. This is not true for the dense phase conveying, where 
the particle interaction is important and where particle density waves can 
be observed. Two modes of dense phase conveying can be distinguished, de- 
pending on the orientation of the pipe transporting the grains. In the hori- 
zontal transport the granular medium separates into two layers, the moving 
bed at the bottom of the pipe and the irregular transport by traveling dunes 
or slugs [9,10,11,12,13]. In the vertical transport the granular medium forms 
strands at medium gas velocities [14,15,16,17,18] and plugs at low gas ve- 
locities with high mass loads [19,20,12]; this mode is called plug conveying. 
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Currently plug conveying is gaining importance in industry, because it causes 
a lower product degradation and pipeline erosion than dilute phase conveying. 

Unfortunately, current models [21,22] of plug conveying disagree even on the 
prediction of such basic quantities as the pressure drop and the total mass 
flow, and these quantities have a great impact in the industrial application. 
One of the reasons for the lack of vahd models is that it is difficult to study 
plugs experimentally in a detailed way. Usually experimental setups are lim- 
ited to the measurement of the local pressure drop, the total mass flux and the 
velocity of plugs. Most promising experimental studies have been performed 
by electrical capacity tomography [20,12] and stress detectors [23,24,13]. Sim- 
ulational studies are handicapped by the high computational cost for solving 
the gas flow and the particle-particle interaction, and are therefore mostly 
limited to two dimensions. For the dense phase regime simulations have been 
done for bubbling fluidized beds [25,26,27,28,29,30], which show for high gas 
velocities first signs of pneumatic transport [31,32,33], for the strand type of 
conveying [16,15,18,34], and slugs in horizontal transport [9,35]. Simulations of 
vertical transport which show plugs with considerable length are rare. There 
is only one simulation of a short pipe in two dimensions by Tsuji et al. [19] 
and a picture of a single plug from a simulation with 50 particles in three 
dimensions by Ichiki et al. [36]. 

The goal of this paper is to provide a detailed view of plugs, by using a 
discrete element simulation combined with a solver for the pressure drop. 
This approach provides access to important parameters like the porosity and 
velocity of the granulate and the shear stress on the wall at relatively low 
computational cost. Contrary to the experiments, it is possible to access these 
parameters at high spatial resolution and without influencing the process of 
transportation at all. Additionally to plug profiles, characteristic curves of the 
pressure drop and the influence of simulation parameters are measured. The 
simulation results are compared with and verified by experiments. 



2 Simulation Model 



Plug conveying is a special case of the two phase flow of grains and gas. It 
is therefore necessary to calculate the motion of both phases, as well as the 
interaction between them. In the following, we explain how our algorithm 
treats each of these problems. 
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2.1 Gas Algorithm 



The model for the gas simulation was first introduced by McNamara and 
Flekk0y [37] and has been implemented for the two dimensional case to sim- 
ulate the rising of bubbles within a fiuidized bed. For the simulation of plug 
conveying we developed a three dimensional version of this algorithm. 

The algorithm is based on the mass conservation of the gas and the granular 
medium. Conservation of grains implies that the density pp of the granular 
medium obeys 

^ + V-(%) = 0, pp = ps{l-<P), (1) 

where the specific density of the particle material is ps, the porosity of the 
medium is (i.e. the fraction of the space available to the gas), and the 
velocity of the granulate is u. 



The mass conservation equation for the gas is 

+ V • (VgPg) =0, OC 0P, (2) 



9pg 



dt 

where pg is the density of the gas and Vg its velocity. This equation can be 
transformed into a differential equation for the gas pressure P using the ideal 
gas equation, together with the assumption of uniform temperature. 

The velocity Vg of the gas is related to the granulate velocity u through the 
d'Arcy relation: 

-VP = ^0(t/,-«), (3) 

where 77 is the dynamic viscosity of the air and k, is the permeability of the 
granular medium. This relation was first given by d'Arcy in 1856 [38]. For the 
permeability k the Carman-Kozeny relation [39] was chosen, which provides a 
relation between the porosity 0, the particle diameter d and the permeability 
of a granular medium of monodisperse spheres. 

After linearizing around the normal atmospheric pressure Pq the resulting 
differential equation only depends on the relative pressure P' (P = Pq + P'), 
the porosity (p and the granular velocity u, which can by derived from the 
particle simulation, and some constants like the viscosity 77: 

dP' Po^, Po 



V{K{(f))WP') - -^Vu. (5) 



dt ri(p 
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This differential equation can be interpreted as a diffusion equation with a 
diffusion constant D = (f)K,{(j))/ri. The equation is solved numerically, using a 
Crank-Nickelson approach for the discretization. Each dimension is integrated 
separately. 

The boundary conditions are imposed by adding a term =pS' on the right hand 
side of equation (5) at the top and the bottom of the tube, where S oc VgPo. 
This resembles a constant gas flux with velocity Vg at a pressure Pq into and 
out of the tube. 

2.2 Granulate Algorithm 

The model for the granular medium simulates each grain individually us- 
ing a discrete element simulation (DES). For the implementation of the dis- 
crete element simulation we used a version of the molecular dynamics method 
described by Cundall [40]. The particles are approximated as monodisperse 
spheres, rotations in three dimensions are taken into account. 

The equation of motion for an individual particle is 

VP 

mx = mg + Fc j- -r, (6) 

where m is the mass of a particle, g the gravitation constant and Fc the sum 
over all contact forces. The last term, the drag force, is assumed to be a volume 

— * 

force given by the pressure drop VP and the local mass density of the granular 
medium — 0). 

The interaction between two particles in contact is given by two force com- 
ponents: a normal and a tangential component with respect to the particle 
surface. The normal force is the sum of a repulsive elastic force (Hooke's law) 
and a viscous damping. The tangential force is proportional to the normal 
force (sliding Coulomb friction) or a viscous damping. The viscous damping 
is used when the relative surface velocity of the particles in contact is small. 
The same force laws are considered for the interaction between particles and 
the tube wall. 

2.3 Gas- Grain Interaction 

The simulation method uses both a continuum and a discrete element ap- 
proach. While the gas algorithm uses fields, which are discretized on a cubic 
grid, the granulate algorithm describes particles in a continuum. A mapping 
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is needed for the algorithms to interact. For the mapping a tent function F{x) 
is used: 

F{x) = mmm, m = - '^/^i' i^/^' - ^' (7) 

[0, 1 < \x/l\, 

where / is the grid constant used for the discretization of the gas simulation. 

For the gas algorithm the porosity 0j and the granular velocity Uj must be 
derived from the particle positions x,i and velocities f j, where i is the index of 
particle and j is the index for the grid node. The tent function distributes the 
particle properties around the particle position smoothly on the grid: 

(t>j^l-Y. Pi^i - ^i)' = rrX H ViF{xi - Xj), (8) 

where Xj is the position of the grid point and the sum is taken over all particles. 

For the computation of the drag force on a particle the pressure drop VPj and 
the porosity (pi at the position of the particle are needed. These can be obtained 

— * 

by a linear interpolation of the fields VPj and 4>j from the gas algorithm: 

0i = ^ (i>jF{xj - Xi) , ^Pi^Yl ^PjFi^J -Si), (9) 

3 3 

where the sum is taken over all grid points. 
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3 Experimental Results 
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Fig. 1. Sketch of the experimental setup. Gas is inserted at a constant mass rate 
at the bottom of the tube, the granulate enters from the left. SI and S2 denotes 
pressure sensors. A CCD camera is used to record the transport at mid-height. 

In the experiment we study the vertical pneumatic transport of wax beads of 
diameter d — 1.41mm, density ps — 937kg /m^ and a Coulomb coefficient of 
0.21. The experimental transport channel is a vertical tube (PMMA) of length 
/ = 1.01m and of internal diameter Df = 7mm. These parameters were chosen 
to be able to compare to simulations with little computational effort possible. 
The air is injected at a constant mass flow rate at the bottom of the tube. The 
granular medium is injected from the side from a tube with a slope of 45°. 

The wax beads, when leaving the vertical tube at the top, are gathered and 
weighted. Prom this data the mass flow rate of the granular medium can be 
obtained. The total pressure drop is detected through two pressure sensors, 
one at the bottom (SI) and one at the top (S2). The velocity and length of the 
plugs can be obtained by a CCD-camera (30 frames per second) positioned at 
the middle of the tube. 
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Experimental results for vertical plug conveying in a tube of length Z = 1.01 m and 
diameter Df = 7 mm using wax beads {d = 1.4 mm) as granular medium. 

Plug conveying is observed, as shown in figure 2. The range of parameters 
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Fig. 2. A series of photos showing a plug moving upwards. The height of the shown 
tube is 9 cm, the frame rate is 30 Hz, from measurement series B. 

for plug conveying in the experiment were rather Umited. The plugs often 
get stuck in the tube, interrupting the plug transport completely. For long 
runs and repeated use of grains, electrostatic forces influenced the results. 
Measurement series for three different air flow rates were made. The results 
for the total pressure drop and the granular flow are given in table 1, for more 
details see [41]. 
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4 Simulational Results 



The setup for the simulation is almost the same as for the experiment. The 
main difference is: particles are injected from the bottom at a constant mass 
flow rate. As default the mass flow rate of the granular medium is chosen 
equal to the measured flow rate in the experiment {2A9 kg/h). Default pa- 
rameters for the particles are: diameter d = 1.4mm, density ps = 937 kg/m^, 
Coulomb coefficient n — 0.5 and restitution coefficient e = 0.5. Simula- 
tions arc preformed for half tube length = 0.525 m and tube diameter 
Df = 7 mm. The resulting gas volume has been discretized into 150x2x2 grid 
nodes, which corresponds to a grid constant of 3.5 cm. The gas pressure is 
set to Pq = 1013.25 /iPa, Simulations are preformed for gas viscosities r] from 
0.045 cP to 0.085 cP and gas flows V between 1.1 l/min and 9.2l/min. The gas 
flow is usually given as the superficial gas velocity Vg — (f)Vg — AV/ttD^ [42], 
which is the equivalent gas velocity for an empty tube. 




time 



Fig. 3. Spatio-temporal image of the porosity along the vertical tube. The dark 
regions correspond to low porosity, light regions to high porosity. The entire tube 
(52.5 cm) is displayed, the elapsed time is 4 s. Default parameters for the simulation 
were used, the superficial gas velocity is Im/s, the gas viscosity is 0.0673 cP, the 
Coulomb coefficient is 0.5. Particles are introduced at the bottom and on average 
3600 particles are within the tube at a pressure drop of 39hPa/m. 

The flow in the experiment is turbulent (particle Re ^ 65). In the simulation, 
an effective gas viscosity is used to account for the turbulence. For an effective 
gas viscosity rj = 0.0673 cP, a gas flow V = 2.31 /min and a Coulomb coeffi- 
cient /i = 0.5 plug conveying is observed as shown in flgure 3. The observed 
pressure drop 39 hPa/m, the plug velocity 0.24 m/s and the plug length 2-4 cm 
flt well with the experimental results. 

As pointed out in the introduction, only very limited simulations of vertical 
plug conveying have been published. Vertical transport has been studied by 
Tsuji [19], but his simulation was limited to a short pipe in two dimensions 
{Lt ~ 3d). In his case the tube contained about 300 particles, which is in our 
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case about the number of particles within one phig. Three dimensional plugs 
of about 50 particles have also been found by Ichiki [36], but no analysis of 
them was done. Thus our simulations are the first to show full featured vertical 
plug conveying at the level of grains in three dimensions. 



4.1 Characteristic curves 



The "characteristic curves" of a pneumatic transport system are plots of the 
pressure drop against the superficial gas velocity Vg = fpVg for different mass 
flows of the granulate. This kind of diagram is highly dependent on the ma- 
terial characteristics of the tube wall and the granulate and can be used to 
predict the overall transport performance for given parameter sets. Such a 
diagram, from data from our simulation, is shown in figure 4. 
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Fig. 4. Total pressure drop against superficial gas velocity for different granular 
mass flows. Plotted are characteristic curves from the simulation for the granular 
mass flows 2.49 kg/h and 9.95 kg/h, and the three data points from the experiment 
with granular mass flows 1.7 — 2.9 kg/h. For mass flow 2.49 kg/h plugs stick when 
the superficial gas velocity is below 0.75 m/s. 

The diagram provides the typical qualitative behavior for pneumatic trans- 
port. Prom top left to bottom right with increasing superficial gas veloc- 
ity, four regions can be distinguished. First, for small superficial gas veloc- 
ities {vs < 0.5 m/s), there is bulk transport. The tube is completely filled with 
granulate, so the pressure drop is high. Nevertheless the drag force on the bulk 
is too small to support the weight of the granular medium. In this case the 
transport comes through the enforced granular mass flow at the bottom. 

For slightly higher gas velocities (0.5 m/s < Vg < 0.75 m/s) the first plugs 
appear. These plugs are often not able to compensate the force from material 
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Fig. 5. Spatio-temporal image of the porosity along the vertical tube. The dark re- 
gions correspond to low porosity, light regions to high porosity. The entire tube 
(52.5 cm) is displayed, the elapsed time is 4 s. A lower superficial gas velocity 
(0.75 m/s) has been used than in figure 3. On average 4900 particles are within 
the tube at a pressure drop of 62 /iPa/m. Due to the low drag force the plugs 
quickly coalesce into very long plugs. 

falling onto them. They usually collapse into big plugs sticking in the tube or 
their motion is reduced considerably (fig. 5). In this region the pressure drop 
decreases rapidly with increasing gas velocity. 

For moderate velocities plug conveying is observed 0.75 m/s < f s < 4 m/s 
(fig. 3). The qualitative behavior of the granular medium in this regime is 
the same in the experiment and the simulation. The particles injected at the 
bottom organize into plugs. After a short acceleration at the bottom plugs 
move upwards with a constant velocity. When two plugs collide, they combine 
into a single plug. A plug always looses particles at the bottom and usually 
maintains its length by collecting particles at the top. A plug disintegrates 
when it gets too small. The porosity of the granular medium within a plug 
is close to the minimum porosity, the edges are sharp and the space between 
plugs is rather empty {(f) > 90%). As one can sec in figure 4 the experimental 
data were measured for a superficial gas velocity at the lower limit for plug 
conveying. In this region the pressure drop is highly dependent on the material 
parameters {ri,Po,fJ^)- 

For high superficial gas velocities {vg > ^m/s) the tube is almost empty 
(fig. 6), in the simulation the particles are pushed out as small plugs. Their 
porosity increases with the superficial gas velocity. In this region the simulation 
method underestimates the pressure drop, because it does not consider the 
increasing drag force on single particles, which in the experiment dominates in 
this region. The regions as described above shift with changes of the simulation 
parameters. 

A nearly proportional relation is observed between the total pressure drop 
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Fig. 6. Spatio-temporal image of the porosity along the vertical tube. The dark 
regions correspond to low porosity, light regions to high porosity. The entire tube 
(52.5 cm) is displayed, the elapsed time is 4 s. A higher superficial gas velocity 
(4m/s) has been used compared to the spatio-temporal image in figure 3 (Im/s). 
On average 850 particles are within the tube at a pressure drop of 13hPa/m. 

and the total number of particles in the tube. This can be explained by the 
observation that the amount of particles between the plugs is small and most 
particles are densely packed at a well defined porosity in the plugs. Through 
the d'Arcy's law, the total pressure drop depends hnearly on the tube length 
filled with this porosity. The pressure drop on the granulate between the plugs 
is negligible and causes only a small deviation from the proportional relation. 
At high gas velocities this is no longer true, because there is no dominating 
porosity for the plugs. 

In the following the dependence of the pressure drop on the air viscosity f], the 
atmospheric pressure Pq and the Coulomb coefficient fi is discussed. For the 
parameter studies the superficial gas velocity has been fixed to Im/s. This 
resembles the gas velocities used for the experiment. For higher velocities the 
parameter dependence decreases. 

The gas algorithm can be influenced by changing the atmospheric pressure Pq 
or the diffusion constant D. An increase in background pressure Pq combined 
with a proportional increase in superficial velocity leaves the pressure drop 
unchanged. This can be deduced directly from equation (5) by noting that 
both terms on the right hand side are proportional to Pq. 

The diffusion constant D oc d^/t] can be changed through the particle diameter 
d and the viscosity rj. Therefore it is sufficient to analyze the parameter space 
for the viscosity at a constant diameter as shown in figure 7. 

The viscosity dependence of the pressure drop can be described by a hyperbola. 
The pressure drop converges for high viscosities to a finite value. In this case 
the pressure drop supports the weight of the granulate embedded in the fluid 
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Fig. 7. Dependence of the pressure drop on the dynamic viscosity r] of the gas at 
superficial gas velocity Im/s. The data points can be fitted by a hyperbola. 

and depends only on the velocity of the gas and the enforced mass flow of the 
granular medium. For small viscosities plug conveying starts to break down 
to bulk transport. The pressure drop diverges at the viscosity where the drag 
force can no longer balance the weight and friction of the bulk. 

The parameter of the particle simulation with most influence on the transport 
of the granular medium is the Coulomb coefficient, i.e. an increase of the resti- 
tTition coefficient from e = 0.5 to e = 0.99 leaves the pressure drop unchanged. 
As one can see in figure 8 the pressure drop increases with raising Coulomb 

45 





0.2 0.3 0.4 
Coulomb coefficient 



0.6 



Fig. 8. Dependence of the pressure drop on the Coulomb coefficient, the superficial 
gas velocity is 1 m/s. For high Coulomb coefficient {fj, > 0.5) a transition from plug 
conveying to bulk transport can be observed. 

coefficients. For low coefficients, < /i < 0.5, plug conveying is observed. 
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for higher Coulomb coefficients the plugs get sticky and coalesces into a big, 
slowly moving plug. These plugs are similar to the ones for low superficial gas 
velocities (fig. 5). 

4-2 Plug statistics 

The spatio-temporal image of the porosity along the vertical tube (fig. 3) pro- 
vides a rough picture of plugs and their movement along the tube. A statistical 
approach is necessary to get more precise values. Properties of interest are the 
porosity and the granular velocity within a plug, the plug length, and their 
dependence on the vertical position z of the plug within the tube. To get 
some average values for the porosity and the granular velocity, the tube was 
segmented into horizontal slices of height 3.5 mm. For each slice the average 
porosity and granular velocity has been computed every 0.01 s. The contribu- 
tion of a particle has been weighed by the volume occupied by that particle 
within a given slice. 

The resulting vertical porosity has been used to identify plugs. Every region 
with a porosity lower than 0.6 is defined to belong to a plug. 




2 I 1 : 1 1 lJ 

0.1 0.2 0.3 0.4 0.5 

z(m) 

Fig. 9. Number of plugs per time as a function of height z. The corresponding 
spatio-temporal image is shown in figure 3. The total tube height is 52.5 cm, default 
parameters are used. The data is averaged over 23 s. 

Figure 9 shows the number of plugs per time as a function of height. At the 
bottom of the tube the incoming granular medium forms into a lot of small 
plugs. Even though their velocity is low, the resulting plug rate is high. As 
can be seen in figure 3, plugs collide along the tube and merge into bigger 
plugs. The merging of plugs reduces the plug rate along the tube. The smooth 
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decrease of the plug rate indicates that there is no preferred height for the 
coUision and merging of plugs. 

For each plug the center of mass, the minimal porosity, the maximal granular 
velocity and the plug length have been computed. 

0.56 
0.54 

0.52 

-e- 

I 0.5 
a, 

0.48 
0.46 
0.44 

0.1 0.2 0.3 0.4 0.5 

z(m) 

Fig. 10. Minimal porosity found in plugs at height z, corresponding to the plug rate 
shown in figure 9. The bars denote the uncertainty of the average, and the dotted 
lines indicate the width of the distribution; at each value of z, half of the observed 
plugs have a porosity lying between the upper and the lower dotted lines. The plug 
porosity decreases at the bottom when the granulate enters the tube (z < 0.1m) 
and increases at the top of the tube {z > 0.4 m), when the preceding plug has been 
removed. 

Figure 10 shows the minimum plug porosity as a function of the vertical 
position 2; of a plug. At each value of z, the mean porosity and its uncertainty 
(standard deviation divided by the square root of the number of plugs) were 
calculated. These quantities are shown by the bars in Figure 10. To show the 
distribution of porosity about the mean, the two dotted lines were added. At 
each height z, half of the plugs have a porosity lying between these two lines. 
The same analysis was carried out for the data in Figures 11 and 12. As one can 
see in figure 10 at the left of the graph, the granular medium is inserted at the 
bottom of the tube with a porosity of 0.51. From there the porosity decreases 
quickly to about 0.47 at a height z = 0.1m and then remains almost constant 
until z = 0.3 m. At the top of the tube {z > 0.3 m) the porosity increases until 
the grains leave the simulation space at z = 0.525 m. 

Figure 11 shows the corresponding particle velocity. The change in porosity 
comes along with an increase of the granular velocity within the plugs. The 
granulate is inserted at the bottom with an initial velocity of 0.04 m/s. The 
granular velocity saturates to a final velocity for the granular velocity at a 
height of 0.1 m. At a height of about 0.35 m the granular medium accelerates 
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Fig. 11. Maximal granular velocity found in plugs at height z, corresponding to the 
plug rate shown in figure 9. The bars and dotted lines have the same meaning as 
in figure 10. The granular velocity saturates in the middle of the tube, the increase 
of the granular velocity at the bottom and at the top of the tube is due to the 
boundary effects described in figure 10. 

until the grains leave the tube. 

An explanation for the final constant plug velocity can be derived using the 
balance equation for the forces on a plug: 

F=-mg-F, + a{(l)){vg-u). (10) 

where m is the mass of a plug, g is the gravity constant, Fc is the force on the 
plug through the friction with the wall and the collisions of particles falling 
onto the plug. The last term a{(f)){vg — u) is the drag force on the plug, which is 
proportional to the relative velocity between the particles and the gas within 
the plug. At the middle of the tube the porosity of the plug (f) is constant. 
On small time scales the mass m and the gas velocity Vg can be assumed 
constant. Only the drag force depends on the granular velocity u through 
d'Arcy's law. For a certain granular velocity the drag force balances the sum 
of gravitational and friction forces. If there is no acceleration of the plug, then 
u remains constant. The solution is stable under small fluctuations in u. 

Figure 12 shows the mean plug length along the tube. The average plug length 
is increasing along the tube, due to merging of plugs. Unlike the previous cases 
no equilibrium length is found. The increase of the plug length combined with 
the decrease of the number of plugs per time (fig. 9) conserves the mass flux 
of the granulate. 

As can be seen in figure 3, the change in porosity, velocity and plug length at 
the top of the tube can be explained by the lack of grains faUing onto the top 
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Fig. 12. Mean plug length at tube height z, corresponding to the plug rate shown 
in figure 9. The bars and dotted lines have the same meaning as in figure 10. The 
average plug length increases along the tube, consistent with the decrease in the 
plug rate (fig. 9). 

of the uppermost plug. 

The diagrams in figures 10 to 12 imply, that there is a typical porosity, granular 
velocity and length of plugs for a given height and a characteristic plug profile 
exists. In the following averaged vertical and radial profiles of plugs at the 
height of 0.26m are discussed. To get some sensible profiles the plugs were 
selected with length and granular velocity close to the mean values (L^ = 
0.024 ± 0.01m, u = 0.16 ± 0.02 m/s). 
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4-3 Vertical plug profiles 



While in experiments recognizing and measuring parameters for global plug 
conveying are rather simple, the measurement of profiles for individual plugs 

remains a nearly impossible task. So one of the reasons for simulating plug 
conveying is to provide a detailed picture of what happens within a plug, and 
how parameters like porosity, granular velocity or shear stress change along 
the plug. 

A porosity profile of an averaged plug is shown in Figure 13. The provided 
profiles (fig. 13-19) average over seven plugs. These plugs were taken from the 
middle of the tube z — 0.26 m with granular velocity 0.16 ± 0.02 m/s and 
plug length 0.024 ± 0.01 m. The coordinate Az denotes the relative vertical 
position along the tube with respect to the center of mass. At the top of the 
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Fig. 13. Vertical porosity profile along an averaged plug, containing about 320 par- 
ticles, positioned at the middle of the tube. The profile was averaged over seven 
plugs with granular velocity 0.16 ± 0.02 m/s and plug length 0.024 ± 0.01m. The 
data belong to the simulation displayed in figure 3. The horizontal axis denotes the 
relative vertical position Az along the tube with respect to the center of mass. The 
vertical position is given in multiples of the particle diameter d = 1.4 mm. 

plug, on the right hand side of figure 13, the porosity decreases quickly from 
90% to 45%. In the middle the porosity of the plug remains almost constant, 
and increases slowly at the bottom. The minimum porosity is bigger than a 
random dense packing due to near wall effects. 

Figure 14 displays the velocity profile of the granular medium in the region 
of the plug shown in Figure 13. One can distinguish four different regions: 
Above the plug (15 < Az/d) the granulate particles are falling downwards, 
the porosity is high. These particles originate from the bottom of an upper 
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Fig. 14. Velocity of the granular medium along the averaged plug of figure 13. 
Inside the plug {\Az/d\ > 7) the velocity is constant, outside the granular medium 
accelerates downwards. 

plug. The drag force cannot support their weight due to the high porosity. 
At the top of the plug the particles collide with the low porosity region and 
decelerate rapidly (7 < Az/d < 15). The granular medium gets denser as 
particles are added to the plug. In the following, this region at the top of the 
plug will be called collision region. Inside the plug, where the porosity settles 
to a low value, the granulate velocity is almost constant {\Az/d\ > 7). At 
the bottom (—10 < Az/d < —7) of the plug the granular medium starts to 
accelerate downwards again. Thus the plug is always loosing material at the 
bottom. This region will be called disintegration region. 

The trajectory of a single particle through the plug can be sketched by a 
snapshot of the vertical forces within a plug. Figure 15 displays the drag force 
Fd and the sum over the interparticle, gravity and friction forces on a particle 
(here called particle forces) Fp. These forces are averaged over the horizontal 
plane. Above the plug, the particle accelerates constantly due to gravity. The 
drag force on a single particle can be neglected, as the particle velocity is 
far from the terminal velocity. The high peak in the particle force marks the 
collision region. The huge fluctuations in the particle force arise, when the 
falling particle collides with the top of the plug front. Within the collision 
region the porosity decreases and therefore the drag force increases. Below the 
collision region, within the plug, the drag force and the particle force balance 
each other. At the bottom of the plug the friction force drops to zero and 
gravity and drag force are in equilibrium. Below this point the disintegration 
region begins, the granular velocity within the plug decreases, the porosity 
increases. After a short distance the drag force is zero and the particles are 
falling freely again. 
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Fig. 15. Forces acting on particles along the averaged plug in figure 13. Interparticle 
forces, gravity and friction are summed up and displayed as particle force. The drag 
force corresponds to the pressure drop of the gas. Outside of the plug (lAz/dl > 15) 
the particle force is mg, indicating that the particle arc freely falling. The peak in 
the particle force corresponds to the collision region, where particles fall onto the 
plug. 




Fig. 16. Stress between the wall and the granular medium along the averaged plug 
in figure 13. The wall stress is limited to the plug, where the porosity is low (< 70%). 

The normal wall stress corresponding to the plug in Figure 13 is shown in 
Figure 16. Due to the implementation of the simulation and the high damping 
the normal wall stress is proportional to the vertical shear stress (Coulomb 
friction). The wall stress is limited to the plug including the collision and 
disintegration region. The wall stress increases from zero within the colli- 
sion region and then slowly decreases, reaching zero in the disintegration 
region. The maximal wall stress is of the same order of magnitude as the 
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pressure imposed on the plug front by the falhng granulate P — ps{l — (f))v^ — 
937 kg/m^ ■ 0.07 • {0.7 m/sf f« 32Pa. 
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Fig. 17. Granular temperature along the averaged plug in figure 13. The granular 
temperature is only non-zero near the collision region between the plug and the 
falling particles at the top of the plug. Due to the large damping (e = 0.5) the 
granular temperature is dissipated quickly and remains on a low level within the 
plug {\Az/d\ > 7). 

Figure 17 shows the granular temperature along a plug. The granular tem- 
perature is the average kinetic energy of particles minus the kinetic energy of 
the motion of their center of mass. Appreciable granular temperatures occur 
in the collision region. Due to the high damping (e = 0.5) these temperatures 
are confined to a small region at the top front of the plug. There is a finite 
granular temperature within the plug, the temperature drops nearly to zero 
at the bottom of the plug. 

4-4 Radial plug profiles 

Another point of view on the plug is given by the radial profiles. The volume 
used for averaging at a given radius is a cylindrical ring with a height of 
1.5 particle diameters and a width of 0.25 particle diameters. The porosity 
at a given range is defined to be the volume fraction of the gas within the 
cylindrical ring. 

As a result of this high radial resolution it is possible to resolve layers of 
granular particles. As can be seen in figure 18, the plugs, averaged for the 
vertical profiles, are highly ordered in the radial direction. The radial profiles 
are taken from the middle of the averaged plugs, at half height of the tube 
{z — 0.26m). The porosity minima in figure 18 correspond to three layers 
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Fig. 18. Radial porosity profile through an averaged plug at the height z = 0.26 m. 
The same plugs were used for averaging as for the vertical profiles in figure 13. 
The radial distance r from the middle of the tube is given in units of the particle 
diameter d = 1.4 mm. The porosity minima correspond to three layers of particles 
moving upwards, parallel to the tube wall. 

of particles. This agrees well with the ratio between the tube and particle 
diameter, which implies that five particles fit into the tube horizontally. The 
minimum at r/d = corresponds to a vertical string of particles in the middle 
of the tube, the layer at r/d = 2 touches the wall of the tube. For higher radii 
r the porosity increases to one. The large error at the center of the tube is 
due to the smaller volume involved for averaging. 
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Fig. 19. Radial profile of the granular velocity for an averaged plug in figure 18. The 
profile shows three layers of particles moving upwards with different velocities. The 
layer along the wall is the slowest due to the friction with the wall. 

The radial profile of the granular velocity (fig. 19) shows that the layers ob- 
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served in the porosity profile move with different velocities. The layer at the 
wall, which has the largest impact on the total averaged granular velocity, 
moves slowest, due to the friction with the wall. 
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5 Conclusion 



In this paper a simple model [37] with coupled grain and gas flow has been 
applied to pneumatic transport. The implementation is three-dimensional; ro- 
tation and Coulomb friction are taken into account. The flux of gas and grains 
is set by the boundary conditions. Plug conveying is observed. The simulation 
used for plug statistics and proflles contained on average 3600 particles in a 
tube of length 0.525 m and diameter 7 mm. The simulation lasted 23 s; during 
this time 90 different plugs were observed at the middle of tube, each with 
about 300 particles. Additionally 24 simulations were preformed to obtain the 
characteristic curves. 

The behavior of the granular medium and the gas is consistent with experi- 
mental observations. The pressure drop, mass flux and superflcial velocity can 
be obtained in quantitative agreement with experiment by introducing an ef- 
fective viscosity rj and an effective friction yU. The effective viscosity reflects the 
increased momentum transport in the gas due to the turbulent flow around 
the grains. The effective friction reflects the complex interplay between sliding, 
rolling and static friction. By using the model, large numbers of plugs could 
be studied, and their porosity, velocity, and size were measured as functions 
of height. 

The simulation results imply that the formation of the plugs at the bottom 
of the tube occurs spontaneously. There is a well deflned preferred velocity 
of the plugs, which is independent of the plug size and the tube length. The 
average plug length increases along the tube due to the merging of plugs. To 
conflrm experimentally the constant velocity and the increase of plug length, 
statistics of plug length and velocity at different heights along the tube are 
needed. The results also show that boundary effects at the bottom and the top 
of the tube have to be taken into account. For the experimental parameters the 
influence of the boundaries is limited to a distance of 0.1 m, which is about the 
distance between two plugs in the middle of the tube. Plugs accelerate in these 
boundary regions. The acceleration at the upper boundary arises due to the 
lack of grains falling on the uppermost plug. This implies that the momentum 
transfer due to falling grains can not be neglected. A model of plug conveying 
must take into account that these falling grains reduce the plug velocity. 

Also a detailed view into an average plug was presented. In contrast to ex- 
perimental setups were are not limited to measure the plug properties only 
at few locations. Experimental results usually provide the plug properties as 
a function of time at a given position along the tube, with the disadvantage 
that the plug proflle is distorted by the relative motion of granulate along the 
plug. The vertical proflles were given for porosity, granular velocity, interpar- 
ticle and drag forces, wall stress and granular temperature. In the experiment. 
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these parameters usually are not accessible along the whole tube. 

Additionally the radial profile is given at high resolution. The results show 
that the grains are ordered in radial direction. The grains arrange themselves 
in layers along the tube wall. These layers have different velocities (up to 40% 

difference). Experimental measurements of the granular velocity using a CCD- 
Camera must take into account that they measure the velocity of the outer, 
slowest layer. 

In conclusion, the model presented here is a useful tool for investigating plug 
conveying. It is fast and flexible enough to make parameter studies on full 
featured plug conveying and provides at the same time access to the plug 
properties at the level of grains. It can obtain both characteristic curves and 
plug profiles. Note also that the model can be easily used to describe hor- 
izontal transport, simply by rotating the gravitational acceleration. Future 
optimizations and faster computers will permit this model to be applied to 
industrial-sized systems. 
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